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TECHNICAL MEMORANDUM X-64910 


SPLASH PROGRAM EVALUATION OF SRB DESIGN 
BACKGROUND 


The Shuttle Solid Rocket Booster (SRB) is designed for recovery at sea 
and reuse. When it impacts the water, it experiences loads that are highly 
dependent upon the conditions of water impact ( velocity and angle) . These con- 
ditions are a function of the meteorological conditions in the impact zone. 
Extremely high winds and seas have a certain probability of occurrence and 
designing for these worst-case loads potentially causes severe cost and weight 
penalties on the structure. Since the vehicle is not manned at time of impact, 
failures are purely economic. Designing for the worst-case load causes high 
program cost due to the resulting conservative design. Designing for a very 
frequent load results in a high program cost because of the large number of 
replacements required with the high attrition rate. Theoretically, the lowest 
program cost results from a lower per unit cost structure than the one designed 
for worst case and a lower attrition rate than that resulting from designing for 
no water impact loads. 

Martin Marietta, under contract NAS8-29622, studied this problem and 
developed a computer program which investigated the attrition rate and costs 
associated with a specific design. The design variables were limited to the 
parachute descent velocity and whether the nozzle extension was jettisoned. At 
that time, the major attrition producing water impact load was slapdown. 

Martin’s computer program was installed on the MSFC computer and 
modified to update the loads and allow for analysis of various designs. It has 
been used to compute attrition rates of seven SRB configurations which are 
described herein. 


PROGRAM DESCRIPTION 


The Martin program is described in Volume II of the report Space 
Shuttle Solid Rocket Booster Recovery Systems Definition. It is a Monte Carlo 
analysis which treats the meteorological i tors (wind, sea, etc.) and the 



strength of each element probabilistically. Each critical load condition is pro- 
grammed as a table of loads input as a function of vertical velocity ( V^) , 

horizontal velocity ( V„) , and water impact angle ( 9 ). For each Monte Carlo 

trial, a water impact condition (V , V , 0) is randomly selected and the set 

V H 

of loads is computed by interpolation from the tables. A set of strengths is 
randomly selected and compared to the loads. If a load exceeds its companion 
strength, a failure is tabulated. The percentage of failures is the attrition rate 
for that particular structure. Motor sinkage is indicated when the slapdown 
load exceeds the strength by a factor of 1. 1 (the ratio of ultimate to yield 
strength for D6AC). 

Using the Martin analysis technique, a new computer program, SPLASH, 
was written which presents the load probabilistically. It presents the entire 
curve instead of the sum of the failures ( the value of the probability) at a 
certain strength level. These results can provide general information about a 
load and they can also be used as a design tool to evaluate a number of designs 
without additional computer runs. SPLASH can present design limit loads 
(predicted actual load with no factor of safety) or it can divide the load by a 
strength ratio probability distribution. The resulting curves are used directly 
with the predicted ultimate strength ( no factor of safety) to determine the 
failure rate. The decision of whether to use the probability of strength is 
based on judgement, Hie failure mode, the k { r<: of testing used, and the degree 
of conservatism desired. 

The meteorological portion of the program has been updated to include 
the latest results of analyses 1 ’ 2,3 . Effects of high altitude wind gusts, low 
altitude wind gusts, wind shear, parachute release dynamics, and wave slope 
are all included. The conditional probability is computed for each contributor 
to the water impact condition so that high gusts will not result with low wind, 
etc. The wave slope computed is filtered to remove the effects of aU waves 
with wave lengths smaller than an effective length unique for each type of load- 
ing. The filtering wave lengths used are vehicle length for slapdown, skirt 
diameter for cavity collapse, and nozzle diameter for nozzle and aft closure 
loads. 


1. SRB Attitude and Horizontal Velocity Probability Distributions at Water 
Impact. S&E-AERO-DD-6-74, MSFC, Mar. 4, 1974. 

2. Natural Environment Inputs for the Monte Carlo Simulation of Sea Surface 
Angle for Shuttle SRB Attrition Studies. S&E-AERO-YA-17-74, MSFC, 

Mar. 27, 1974. 

f . Estimation of SRB Coning Motion for Attrition Studies, S&E-AERO-DD-9-74, 
MSFC, May 3, 1974. 
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The actuator loads are a unique problem because of the azimuth orien- 
tation of the actuator relative to the horizontal velocity vector. To include this 
effect, an azimuth angle was selected randomly and the in-plane actuator load 
was multiplied by the sine of that angle. 

The strength probability 4 is dependent upon the type of testing. Thp 
options for which data are available are no test, standard test, proof test, and 
model test. The results enclosed are for the standard test where a full size 
prototype or structural model is tested to the design limit load and corrections 
to the analytical model and the design are made as a result of any failures. 


RESULTS 


The strength probability distribution for standard test is shown in Figure 
1, and the water impact conditions are shown in Figures 2 and 3. Note that 99 
percent probability values of water impact are that the angle will be less than 
12 deg and the horizontal velocity will be less than 19 m/sec ( 63 ft/sec) . 

Presented in Figures 4 through 1 8 are probability distributions for the 
significant water impact loads. In each case the appropriate wave slope filter- 
ing has been included and the curves with and without probability of strength are 
shown. Cases with and without the nozzle extension are plotted separately. To 
properly use the load curves, enter the curve with a design capability on the left 
and read the attrition rate on the top scale or the probability of nonexceedance 
on the bottom scale. 

Generally, the probability of strength effect is to reduce the probability 
of failure for structures that have been designed for about 90 percent loads or 
less, For structures that have been designed for 3 sigma loads, the probability 
of low strength (less than 10 percent, Figuie 1) will increase the failures when 
it is included and thus the curves tend to cross over on the right side. 

The wave slope is added directly to the angle computed using wind and 
parachute dynamics. The effect of filtering the wave slope on the resulting angle 
can be seen in Figures 3 and 4. The filter length is based on judgement and was 


4. Thomas, Jerrell and Hanagud, S.: Reliability- Based Econometrics of Aero- 
space Structural Systems: Design Criteria and Test Options. NASA TN D-7646, 
June 1974. 
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chosen upon the recommendation of D. Kross of MSFC. It assumes that the wave 
half length must be long enough to contribute to the dynamics and, therefore, 
must approximate the surface length over which the load acts. The effects of 
filtered wave slope on the load distribution are less dramatic than on the angle 
itself. The wave slope effects on all loads except cavity collapse and actuator 
were not visible on the plots. 

A recent development in the meteorological investigations has been the 
introduction of coning effects. The parachute/SRB dynamics are three dimen- 
sional and include a nutational pendulum motion that combines rotation about 
several axes simultaneously. An additional problem is that the coning dynamics 
( rotation about the local vertical axis) have less damping associated with them 
so they increase the probability of higher angles associated with low horizontal 
velocities. Unfortunately, high angle and low velocity is a high load case for 
cavity collapse. Dynamics supplied by the Systems Dynamics Laboratory at 
MSFC 5 have been programmed lor the nozzle and cavity collapse loads. The 
effects of wave slope and coning on the cavity collapse load can be seen in 
Figure 4. The distribution including all effects is plotted as a solid line and the 
load with wave slope and coning effects individually removed are plotted as 
dashec lines. 

The cavitj collapse loading case was a particular problem because the 
pressure distribution changed shape as well as peak \ .lue. Some of the loading 
conditions presented in MSFC document S&E-ASTN-ADL ( 74-38) 6 that had the 
highest peak pressures were not the most severe loading conditions. Because of 
this, the peak pressure (P 2 ) cannot be used as an indicator of failure. Several 
possible indicators were investigated. Lee load (the area under the pressure 
curve on the lee side) shows the best agreement with eigenvalues resulting from 
structural analyses and was used for the attrition rates computed. 

The configurations evaluated are shown in Figures 19 22. Configuration 
0 (Fig. IP) is the proposal configuration. Configuration 1 (Fig. 20) is a design 
that includes capability for cavity collapse. Configuration 1-1, modified (Fig. 

2i ) , has been optimized with the performance margin. Configuration 3 ( Fig. 22) 
is a performance- only configuration that has the performance capability of 


5. See footnote 3. 

6. Updated SRB Cavity Collapse Water Impact Loads, Configuration Without 
Nozzle Extension. S&E-ASTN-ADL ( 74-38) » MSFC, Apr. 26, 1974. 
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configuration 1-J but has no weight or design for water impact loading. Tables 
1 and 2 contain the results of comparing the capabilities of each configuration 
with the probability of loads. The capabilities and attrition rates are tabulated. 
In Table 2, maximum and minimum values are given based upon the following 
assumptions: 

Actuator yield strength is 1112 kN (250 000 lb). 

Actuator ultimate strength is 1334 kN (300 000 lb) (FTU) 7 . 

Nozzle ring yield at 1112 kN (250 000 lb) actuator force. 

At 1334 kN (300 000 lb) , nozzle is released to impact other structure 
resulting in the following attrition rates: 


Throat and Seal 

Skirt 

Closure 

TVC Fluid Loop 


Maximum 

FTU Rate 
0.5 FTU Rate 
FTU Rate 
FTU Rate 


Minimum 

0. 1 FTU Rate 
0.05 FTU Rate 
0. 1 FTU Rate 
0. 5 FTU Rate 


COHERENT DESIGN REQUIREMENTS 


The probability distributions can be used to establish coherent design 
requirements for water impact. Coherent requirements are requirements for 
each load that result in similar failure rates. They may or may not be realistic 
to design to, depending on the cost/attrition trade of each load, interrelation- 
ship between loads, and other factors. The following are coherent design require- 
ments for 1 percent water impact attrition: 


7. FTU means force tension ultimate. 
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Factor of Safety Included 
W/O Extension With Extension 


20 N/cm 2 ( 29 psi) 
17 N/cm 2 (25 psi) 
570 N/cm 2 (827 psi) 
1.33 x 100/15/0 
Distribution 


21 N/cm 2 (30 psi) 
18 N/cm 2 (26. 7 psi) 
46 N/cm 2 (67 ps : l 
1. 65 x 100/15/0 
Distribution 


Slapdown 
Submergence 
Aft Closure Pressure 
Cavity Collapse 

Nozzle Side Load 
Nuzzle Axial Force 


1735 kN (390 000 lb) 
1557 kN (350 000 lb) 


3000 kN ( 675 000 lb) 
11.4 MN (2 570 000 lb) 


ATTRITION PROGRAM 


To use the results of the SPLASH program, another computer program 
was written 8 to compute the total number of units required. This program 
uses the mission model as supplied in the requirements documents. An 
attrition rate function that reduces with time and has an average attrition rate 
as found from the SPLASH program is input. The output includes the number 
of units required for the Shuttle program. The effects of refurbishment time, 
maximum life, and the variable launch rate are all included. 

Tables 3 and 4 contain the total attrition rates and units required for 
445 flights. The total attrition rate is 3 percent (for all nonwater impact 
causes) added to the individual structure failure rate plus the sinkage rate as 
shown in Tables 1 and 2. Where two numbers are given in Table 4, the larger 
number represents the maximum failure assumption for actuator cascading 
failure effects, and the smaller number is the minimum assumption. These 
numbers of required units can then be used in a cost .malysis to compare the 
various designs and to determine the most effective design. 


8. Program was written by the Operations Development Branch of the Systems 
Analysis and Integration Laboratory at WSFC. 
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TABLE 2. ATTRITION RATES FOR ACTUATOR- INDUCED FAILURES 
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TABLE 3. TOTAL ATTRITION RATES (%) 


Configuration 



0 

0-X 

a 

1-X 

liiL 

1.1-X 


Aft Cylinders 

100 

100 

17 

29 


D 


Foward Cylinders 

n 

Bl 

3. (5 


3.4 



Max 

■ 

21 

3.4 

32 

3.3 

36 

SRM 

Aft Closure 







Min 

■ 

5.3 

3.3 

6.3 

3.2 

6.6 


All Other Motor 
Seg. 

3.4 

n 

3.3 

3.4 

3.2 

3.3 


Nozzle Ring 

Qj 

27 

4.5 

41 

4.0 

45 
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32 
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and Seal 
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Figure 1. Probability of strength, standard test. 




HI 99.8 9 M M SO SO 70 60 SO 40 30 20 10 S 2 1 0-: •' 1 9 0S »vi 



Figure 2. Horizontal velocity probability. 
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Figure 3. Water impact angle 












Figure 5. Nozzle side force, extension jettisoned, 
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Figure 6. Nozzle side force, extension retained. 
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Figure 8. Nozzle axial load, extension retained. 
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Figure 9. Aft closure pressure, extension jettisoned. 
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Figure 10. Aft closure pressure, extension retained. 



















Figure 14. Cavity collapse load, extension retained. 











Figure 16 . Submergence pressure, extension retained. 






Figure 17 . Slaptiown pressure, extension jettisoned 









Figure 18. Slapdown pressure, extension retained. 
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Figure 19. Case segment configuration 0. 
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Figure 20. Case segment configuration 1. 
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Figure 21. Case segment configuration 1-1 modified. 





Figure 22. Case segment configuration 3. 
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